The classical J1-J2 Heisenberg model on the Kagome lattice 
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Motivated by earlier simulated annealing studies and materials with large spin on the Kagome 
lattice, we performed large scale parallel tempering simulations on the Kagome lattice for the ex- 
tended classical Heisenberg model including next nearest neighbor interactions. We find that even a 
small inclusion of a J2 term induces anti-ferromagnetic order which prevails in the thermodynamic 
limit. The magnitude of this effect is surprising. While at J2 — the finite-size behaviour does 
not suggest a phase-transition, at other points the numerical result is consistent with one. Close 
to J 2 = and for a positive sign of J 2 two subsequent phase-transitions/crossovers are found, one 
of them connecting to the crossover for the J 2 = case, shedding light to the pure case. The 
universality classes of the transitions were explored. 
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INTRODUCTION 



The 2D Kagome lattice being the underlying lattice 
of a number of materials, is of great interest due to 
its highly frustrated geometry with respect to the 
antiferromagnet. A number of these materials have 
larger spin. E.g. SrCrg-^Ga^^Oig 1 with spin-| 
or KFe 3 (OH) 6 (S0 4 )^ with spin-f. With growing 
spin we expect the model to exhibit some behaviour 
captured in the classical model. This model has been 
studied by Chalker et al 4 by Monte- Carlo, spin-wave 
approximation and low temperature expansion. Here the 
"order-by-disorder" scenario with coplanar ordering was 
proposed^ Reimers et al£ have performed Monte-Carlo 
runs more systematically and could confirm some of the 
predictions, particularly coplanar ordering. Chandra 
et al and Ritchey et al have extended the work and 
studied anisotropy^ A cross-over for the pure case 
and a second-order phase transition for any inclusion 
of anisotropy was proposed. More recently the classical 
Heisenberg model has been investigated by Zhitomirsky^ 
with a simulated annealing method. He found a strange 
finite-size behaviour in the specific heat, inconsistent 
with a phase-transition. The antiferromagnetic order 
was shown to vanish in the thermodynamic limit and the 
so called "octopular-order" parameter, a third rank ten- 
sor, was proposed to describe the low temperature phase. 
Subsequently Gvozdikova et al extended the studies 
by inclusion of a magnetic field in the Hamiltonian 10 
following earlier work of Zhitomir sky j^ 1 - Another relevant 
study is by Korshunov, who studied the XY-Ji — J 2 
model on the Kagome and found magnetically and 
topologically ordered phases* 1 ^ 

Classical models with a continuous symmetry are 
generally not expected to show any long-range order 
at finite temperature, except for topological order, 
due to the theorem of Mermin and Wagner.— For 
frustrated models the situation is more subtle, as a 
discrete symmetry can be induced and this symmetry 



can in principle be broken. This possibility was first 
suggested by P. Chandra et aU^, and has been ob- 
served numerically by Weber et al, 15 . Note that this 
does not violate the Mermin- Wagner theorem, as the 
symmetry being broken is not a continous one, but 
an additional discrete one. Meanwhile a mathematical 
proof is published about such type of models^ showing 
that such induced symmetry can be broken. In the 
classical nearest neighbor Heisenberg-model on the 
Kagome (CKLNNHM), all studies of the past are not 
consistent with a phase-transition at finite temperature. 
Features in the observables, such as a peculiar finite 
size behaviour in the specific heat have been observed 
in Monte- Carlo simulations^ We refer to this point as 
"crossover-temperature" in our manuscript. 

Here we investigate the extended Heisenberg model 
by inclusion of a next-nearest neighbor (NNN) coupling 
term. The motivation for this is three- fold: First, in 
real materials a small NNN coupling is expected; sec- 
ond, from observing the trends towards J 2 — » new in- 
sight to the pure model might be possible; and third, this 
model couples to three underlying sub-lattices which are 
again of Kagome type (see Fig. [T]). This is therefore the 
first study of frustration which considers the coupling 
of two frustrated lattices which to some extent "unfrus- 
trate" them as the inclusion of the J2 lifts some of the 
degeneracies. Our results show that even a very small 
J2-term enables an anti- ferromagnetic order. This order 
is of q = (0, 0) type for AF-coupling of J2, and of \/3x a/3 
type in the other case. For the pure case (J2 =0) the 
crossover temperature is tiny (T& ~ 0.003) this changes 
quite dramatically upon switching on the J2-term (e.g. 
we obtain T c = 0.045 ± 0.001 for J 2 = 0.02). Away 
from J2 = the finite-size scaling is consistent with a 
phase-transition. The universality class could be classi- 
fied to some extent and was found to be of the Ising type 
for J 2 ~ ±0.02. For positive J 2 which is close to J 2 = 
two peaks in the specific heat are found numerically, sug- 
gesting two phase- transit ions. Studying the limiting be- 
haviour from the other side towards J2 = appears to 
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FIG. 1: (Color online) The Kagome lattice, including the 
next nearest neighbor bonds (in dashed lines). The three 
underlying NNN-Kagome sub-lattices are shown in different 
colours. 

be numerically harder. Here the phase- transit ion is con- 
sistent with a first-order transition, which is apparently 
the reason for the numerical difficulty. 

II. MODEL 

From earlier work by the authors mentioned in the in- 
troduction, it was shown that the CKLNNHM is very 
hard for a numerical investigation: The interesting tem- 
perature region is very low, and the low temperature 
phase is very difficult to bring into equilibrium. The 
extension of further neighbor interaction lift some degen- 
eracies and therefore taking away some of the numerical 
"hardness" . Our model Hamiltonian reads: 

H = -J x • Sj + J 2 ^S, • Sj (1) 

where the J\ term couples between nearest neighbor 
pairs (indicated by solid lines in Fig. [Tj) and the J 2 term 
between next-nearest neighbor pairs (indicated by dashed 
lines in Fig. []]). In the following we give J 2 in units of J\. 
S is a classical spin operator. When J 2 = we recover 
the nearest neighbor model. In this case there is a large 
degeneracy of coplanar and non-planar states. Let us 
discuss the possible ground states in this case: Ground 
states can be found by rewriting the Hamiltonian from a 
nearest neighbor sum to a sum over all triangles a as: 

H = -J 1 J2Si-Sj = y^(Si+S 2 +S 3 ) 2 + c (2) 

(ij) a 

where c is a constant term. From this sum one realizes 
the condition to minimize the energy: If the sum of 
each individual triangle is zero, this is a minimum 
energy state. There are two groups of such states: the 
coplanar ones, where all 3 spins lie in one plane, and the 
non-coplanar ones, where this is not the case. We will 
first discuss the coplanar states and briefly mention how 



to obtain the non-planar ones from them. The coplanar 
states contain two types of antiferromagnetic sub-lattice 
order pattens (see Fig. [2] thereafter AFM-states). There 
is a continuous 0(2) symmetry for these states as a 
simultaneous rotation of all spins in the plane keeps the 
energy constant. Both variants of AFM type states are 
two-fold degenerate: For the q = (0, 0) state (from now 
on "AFM I" state) there is one version with only positive 
chiralities and one state with only negative chiralities. 
Chirality refers to the clockwise or anti-clockwise type 
of rotation of spins when circumferencing an elementary 
triangular plaquette. For the other (q = (47r/3,0) from 
now on dubbed as "AFM II" state - this state is usually 
referred to as the y/S x y/3 structure) AFM state there 
is one version with all the up-triangles having positive 
chirality and all the down-triangles having negative 
chirality and one version where the opposite is the 
case. For an illustration please see Fig. [2j The above 
condition for minimizing the energy is not only full-filled 
by the two AFM states, but by all states having 3 spins 
in one plane arranged among the plaquette with 120 
degrees. Starting with one sub-lattice, the angles are 
fixed, but the arrangement of the 3 spins can be varied. 
For an arbitrary arrangement fullfilling the 120 degrees 
condition there is no regular patten in chiralities. One 
may label the 3 spins with 3 colors e.g. red, blue and 
green, and compose pattens where one coloured site has 
only nearest neighbors of another colour. All coplanar 
states can be constructed in such a manner. See Fig. [3] 
for illustration of a general planar state. 

Taking a coplanar states as a starting point, one 
may use deformations leaving the energy unchanged to 
construct non-planar states: If e.g. we consider the two 
AFM states of Fig. [2] there are two types of defects: 
The weather-vane defect for AFM II, here one rotates 
all inner spins of a hexagon along the common axis 
given by the outer spins, and a line-defect for the AFM 
II, where one can rotate a whole line, by the axis given 
by the neighbouring spins. The energy is kept constant 
in both cases. Considering an arbirary situation (a 
coloured state which is neighter AFM I nor AFM II), 
one can do such operations among two coloured lines. 
From the above discussion one realizes that the ground 
state is highly degenerate and disordered. As mentioned 
in the introduction, previous studies proposed an 
order- by-disorder scenario: For T / entropic selection 
favours coplanar states over the other states. As an 
order parameter originally a nematic order parameter 
was considered in the earlier studies. This measures the 
alignment of spins with respect to a common plane. The 
general patten of 120 degree alignment, discussed above 
and being an additional symmetry is not measured by 
such an order parameter. The octopular tensor was 
proposed as a better choice to measure states as the one 
in Fig. 03 We follow this approach. 

Finally we turn to the situation with J 2 7^ 0: In this 





FIG. 3: (Color online) An octupolar ordered state. 



FIG. 2: (Color online) Possible anti-ferromagnetic (AFM) 
sub-lattice ordering pattens. Top: AFM ordering with q = 
(0, 0) with an unit cell similar to the one of the Kagome lat- 
tice, in its two possible forms with positive (top left) and 
negative (top right) chiralities. We refer this state as AFM I 
in the text. Below: AFM ordering with q = (4tt/3,0). The 
corresponding enlarged unit cell is indicated in grey. This 
AFM state possess two possible forms as indicated in the two 
pictures. We refer to this order as AFM II from now on. 



case we bias the lowest energy state depending on the 
sign chosen for J2 towards the two AFM states (which 
are still degenerate with respect to rotation and have 
an additional two- fold degeneracy as discussed above). 
Therefore we lift most of the degeneracies. When being 
in one perfectly ordered AFM state of one manifestation 
of the chirality, one needs to overcome an energy barrier 
to reach the other one. This is different from the situation 
at J 2 = 0, as there the two states can be reached with 
formation of domain walls costing zero energy. There- 
fore at J2 7^ there is a discrete Z2 symmetry for planar 
states. The building blocks of up and down triangles pos- 
sessing and '-' chiralities can therefore be understood 
of taking the role of Ising variables. Note that for planar 
states the situation is the same as in the XY-model on 
the Kagome lattice being discussed in detail in Ref. [TH. 



III. METHOD 

We employ an improved parallel tempering (PT) 
approach^ — based on the ALPS libraries and 
parapack^ We have tested several optimizations 
from the literature and finally used the following 
method: A suitable set of temperatures is chosen, 
then auto-correlation-times are measured to optimize 
the frequency of exchange steps between neighbouring 
replicas. This is the so called "r-method" of Bittner et 



alj22 As updates we used restricted updates and perform 
over-relaxation, 9 This method shows roughly 100 times 
faster convergence compared with canonical Monte 
Carlo and naive updates. Typically we used around the 
order of 10 5 sweeps for thermalization and another 10 7 
sweeps for measurements. One PT simulation consisted 
of 90 to 180 replicas. The exact amount of sweeps for 
thermalization and measurements and the number of 
replica varies largely depending on the regime of interest. 
The AF-coupled region close to J2 = is very hard. 
As discussed below we believe this is a consequence of 
the nature of the phase-transition at this point. We 
have tested our program intensively by comparing our 
data against earlier published studies.— i 22 i 23 The J 2 term 
was tested by noting the structure of the NNN-lattice 
being of 3 independent Kagome lattices of smaller size 
(see Fig [I]). Choosing the couplings J\ = 0, J 2 — 1 
should therefore give similar curves in the measurements 
compared to choosing J\ = 1, J 2 = when considering 
the change in the volume-size N. Our tests left no doubt 
our program works correctly. 

To identify possible phase-transition points, we mea- 
sure the specific heat c v . This is an quantity which should 
show a feature at a transition or even a crossover point, 
independently of the kind of order-parameter: 



(E 2 ) ~ {Ef 
NT 2 



(3) 



Following earlier literature we measure the octupolar 
order parameter Tq P as mentioned above and investi- 
gated by Zhitomrisky&2 



( Ta/37 ) 2 = ^E[(( s - s ,) 3 )-|(^-s,) 



(4) 



with the limiting value \ in the limit T — >• 0. See 
Fig. [3] for an example of an octupolar ordered state. In 
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addition we measure antiferromagnetic long-range order 
defined as: 

where the first index I labels the site within the unit- 
cell, the other indexes label the xy-position of the unit cell 
within the lattice and Ri is the position of the unit-cell 
within the underlying triangular lattice. For q we choose 
the two vectors corresponding to the ordering pattens 
for the AFM I and AFM II, q = (0, 0) and q = (4tt/3, 0). 
Among the two orders the AFM I order has been shown 
to be very difficult to equilibriate for the pure case, par- 
ticular cooling and warming runs starting from this state 
differ significantly. 6 On the other hand there have been 
speculation about possible AFM II order or correlations 
in the pure model^i which have been shown to van- 
ish for any non-zero temperature performed numerically^ 
Another quantity of possible interest is the spin-stiffness: 

ps = -^{\ {E) + y( E (Si x Sj) • ^ • elj) ] 2 )} 

<ij> 

(6) 

where e a and e# are both arbitrary unit vectors, eij 
is the unit vector in the direction of the bond between 
i and j. To investigate possible further structure within 
the spin-ordering, we measured the spin-structure factor 
S(q) for a few particular values defined as: 

S(q) = ^£(S;-S,)e^-^ (7) 

where Ti labels the real space position, and q refers to 
the momenta. 



IV. RESULTS 

Our main result is summed up in the form of a phase 
diagram in Fig. [H As established below at any J2 7^ 
we see a finite-size behaviour consistent with a finite 
temperature phase transition. Starting from the 
crossover temperature at J 2 = 0, we see that a small 
value of J2 does not only enable a phase transition 
but has a huge impact on the point of transition. 
As expected depending on the sign of J2 we achieve 
ordering of the AFM I or AFM II type. The transition 
has been found to be of second order with 2D Ising 
universality class for | J2I > 0.0075 and of first-order for 
> J 2 >= —0.0075. For tiny but positive J 2 two peaks 
in the specific heat are observed and a classification in 
terms of universality failed. 

The finite-size behaviour is summed up in Fig. [5l Here 
we show the three most relevant observables Tq P , tuaf 
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FIG. 4: (Color online) Phase-diagram of the extended 
Heisenberg-model on the Kagome lattice: Shown are the esti- 
mated transition (for non-zero J 2) and crossover (for J2 = 0) 
points corresponding to antiferromagnetic and octopular or- 
der. AFM I labels the phase with the q = (0, 0) antiferro- 
magnetic order, AFM II labels the phase with the x 
structure - in these two phases the octupolar order is found 
in addition to the antiferromagnetic one, IP labels the region 
between the two peaks of the specific heat for the positive 
J2 -branch close to J 2 — 0. The solid line indicates the transi- 
tion line consistent with first order scaling laws. The dashed 
lines are those transition lines consistent with the 2D Ising 
universality class. For the dotted line the classification failed. 



and c v for four values of J2. These are snapshots for 
the interesting regions in the phase-diagram. For values 
of \J 2 \ > 0.0075 (in the picture we choose J 2 = ±0.02) 
there is a phase transition between the disordered and 
AFM phase. The specific heat develops a peak which 
gets narrower and higher with increasing lattice size. 
The AFM phase is of AFM I type for negative J 2 
and AFM II type for positive J 2- The position of the 
peak appears where the measured AFM-ordering has 
its inflection point. At the same position Tq P appears 
to have its onset. For smaller values of J 2 the nature 
of the transition changes: For e.g. J2 = —0.005 the 
peak of the specific heat grows much faster compared 
to the IJ2I > 0.0075 cases. The order-parameters tuaf 
and Tq P have a more sudden change between the two 
phases. In this part there is very slow convergence when 
simulating. To us it appears to be a first-order transition 
which we elaborate below. In the region close to J2 = 
but with positive sign for J2, we find another interesting 
feature: For e.g. J2 = 0.00125 (see Fig. [5|) we find two 
peaks in the specific heat, the taller one matching to 
the inflection point of vtlaf the second, less dominant 
one matching to the onset of the Tq P and tuaf order. 
When observing the finite-size effect, we notice that the 
bigger peak stays almost stationary at the same position 
with respect to the temperature, when increasing the 
system-size, while the smaller peak moves towards 
smaller temperature with increasing system-size. An 
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extrapolation suggests the two peaks to stay separated 
in the thermodynamic limit. Investigating the region 
between the two peaks we found that both Tq P and 
ttiaf orders are reduced in the thermodynamic limit: 
an extrapolation suggests here the Tq P order clearly to 
vanish, while for the AFM a small value may remain for 
infinite lattice sizes. This "phase" therefore seems to 
resemble some of the behaviour compared to the J2 = 
case, which shows similar features in Tq P and m^. 
Two differences of this intermediate "phase" compared 
to the situation at J 2 = can be seen: First, in the 
J2 = case the vtlaf order is vanishing faster than 
the Tq P , here this trend is opposite. Second, the finte 
size trend for the upper specific heat (the smaller one) 
seems still consistent with one at a phase-transition as 
this peak still seem to grow taller and narrower with 
larger system-size. Since the two peaks are quite close 
to each other, it is hard to be clear about the trend of 
this second smaller peak, resting on top of the slope of 
the bigger peak. This upper transition-line separates at 
around J 2 = 0.002 from the main-phase boundary and 
continues to the point of crossover for the J 2 = case. 
Using the transistion points to extrapolate to the regime 
where simulation is no more possible it one notices: The 
lower transition line for the positive J2 branch between 
the AFM 2 state and the disordered phase tend to 
T for J 2 -> 0, while in the other branch the AFM 1 
order seems to vanish before reaching the J 2 = point. 
This is consistent with earlier claims of a T = critical 
point for the AFM 2 transition^ For J 2 = 0.00125 
we checked to which of the two peaks the spin-stiffness 
matches (for the other cases, the spin-stiffness generally 
speaking has the same behaviour in terms of onset and 
infliction point compared to the tuaf and Tq P ). This 
is shown in Fig. [6l We observe again a similar trend 
compared with the Tq P and tuaf measurements: The 
onset matches the first peak, while the infliction point is 
matching the broad peak. 

In Fig. and Fig. [8] we plotted the development of 
the measured quantities for a larger system size (L = 24) 
as a function of J2. For the negative branch, one can 
see clearly that the nature of the transition is becoming 
more consistent with a first-order phase-transition when 
coming closer to the point J 2 = 0. In the positive 
branch, the specific heat is changing from a very broad 
single peak at large J2 values to a smaller double peak 
close to J 2 = 0. The octupolar order Tq P has a more 
sudden rise when tuning J 2 towards smaller values. 
At J2 = the peak in the specific heat vanishes when 
increasing the lattice-size, as observed by Zhitomirsky, 
suggesting no phase transition at this point. The AFM 
order vanishes at this point in the thermodynamic 
limit. Elsewhere, as mentioned above we note that the 
trend of the specific heat, seems to be consistent with a 
phase-transition, this is true also for the second smaller 
peak for positive and small J2. In We tried to fit the 
data to scaling laws: The points easiest to fit, are those 




FIG. 5: (Color online) The finite-size behaviour of the most 
important measurements for four values of J2. Error-bars are 
smaller than the symbol-size. 
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FIG. 6: (Color online) Spin- Stiffness for J 2 = 0.00125 and 
several system sizes. Error-bars are smaller than the symbol- 
size. 



with large values for | J 2 | as here the temperature for the 
point of transition is higher. Fitting a relation for the 
dependence of the estimate of T c on L (we tried several 
ways, one of them is the peak in c v ) to a power law gives 
a value for the critical exponent of v = 0.99 ± 0.05 at 
J 2 = 0.02 where the error-bar is estimated from different 
ways to obtain the T C (L). This would be consistent with 
a 2D Ising universality class. To test this theory we 
fitted the height of the specific heat. For the 2D Ising 
universality class such a fit should not be possible, as 
the divergence is logarithmic. Any power-law fit gives 
indeed systematic errors compared to the data, while a 
logarithmic fit agrees well (see Fig. [9] lower panel). For a 
value of J2 = 0.01 we obtain v = 1.08 ±0.05, and again a 
logarithmic divergence for the height of the specific heat. 
This still seem to somewhat fit in this universality class. 
For smaller values of the positive signed J2, our data 
does not allow to get a consistent value for this critical 
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FIG. 7: (Color online) Measurements for the negative J2 case 
of the antiferromagnetic order (tuaf), the octopular order 
(Tq P ) and the specific heat c v for several values of J2. The 
system size was L — 24. Error-bars are indicated if they exeed 
the symbol size. 
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FIG. 8: (Color online) Measurements for the positive J 2 case 
of the antiferromagnetic order (ttiaf) the octopular order 
(Tqp) and the specific heat c v for several values of J2. The 
system size was L = 24. Error-bars are smaller than the 
symbol size. 

exponent. Here the situation is more complicated since 
two peaks in the specific heat develop. As these two 
peaks are close to each other, it is not too surprising 
that finite size scaling fails. At the negative J2 branch, 
the large | 1 value is again consistent with the Ising 
universality class with a value of v = 0.99 ± 0.05. For 
smaller values of \J<i\ of the same sign, the transition 
point appears to become sharper and the specific heat 
diverges faster, as one would expect in a first-order 
transition. For a first order transition the maximum of 
the specific heat is expected to be proportional to the 
volume L 2 and T C (L) should scale with L~ 2 £± Plotting 
Cv,max versus L 2 for the case J2 = —0.005 gives a linear 
dependence and the scaling of T C (L) agrees well with 
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FIG. 9: (Color online) Scaling fits for two points J2 = —0.005 
(above) and J2 = 0.02 (below). Above left: The specific heat 
height is found to be proportional L 2 ; Above middle: T C (L) 
scales with L~ 2 ; Above right: Double peak structure in the 
energy at the estimated transition point: We used around 
20 mill, sweeps for measurements and the system size was 
L = 24; Below left: T C (L) can be fitted towards power law. 
Below right: Peak height shows logarithmic dependence. 



the expected L~ 2 relation (see Fig. [9]). Therefore we see 
that the maximum of the specific heat is proportional to 
the volume and the position of the same peak is propor- 
tional to L~ 2 as expected for a first order transition. We 
recorded energy-histograms and found the double-peak 
structure a typical feature of first-order transistions. 
Obtaining this histogram was somewhat challenging, 
since it converges very slowly and seem to appear only 
very close to the transition point (as we show in Fig. 
[9]). We needed 6 digits for the temperature to find two 
peaks and 20 mill, sweeps for measurements. From 
our observations we estimate the change from second- 
order to first-order behaviour taking place at around 
J 2 w —0.0075 ±0.002, here a tricritical point is expected. 
There are several possible scenarios how a second-order 
phase transition changes to become a first-order one 
when tuning parameters^ an analysis of the under- 
lying mechanism causing this, we leave for future studies. 

We investigated the magnetic structure factor, to see 
the impact of the inclusion of the J2-term on this quan- 
tity. For large | J2 1 we have chosen two temperatures, one 
above the transition (T = 0.1) and one below (T = 0.02). 
For the small IJ2I cases we chose two more very low 
temperatures T = 0.005 and T = 0.002, the structure 
for these two very low temperatures is the same, there- 
fore we included only T = 0.005 in Fig. [101 As the 
Kagome lattice is not bi-partite, the structure factor ex- 
tents over several Brilouin zones. For high temperatures 
we observe the diffuse Bow-tie structure, as in the pure 
Kagome Lattice, 9 and as in related lattices such as the 
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FIG. 10: Magnetic structure factor for L = 9 systems for several values of J2 and for the relevant temperature range. 



Pyrochlore^ The structure appears for all investigated 
values for J 2 but at different temperatures. Lowering the 
temperature, from there we find qualitatively the same 
structure than for the J 2 = case for positive sign of 
J 2 (3 peaks at the center of the bow-tie prevailing). For 
negative J 2 with the other AFM ordering for the low- 
temperature phase, three peaks at the boundary of the 
triangle of the bow-tie develop, when lowering the tem- 
perature. In both cases the lowest energy plots finally 
show pinch points at the AFM ordering vectors for the 
state concerned but with an increased periodicity, while 
at high temperature the correlations are similar to the 
J 2 = case. This can be seen from our plots in Fig. [TOl 
For e.g. the tiny values of | J 2 | shown (J 2 = —0.000625 
and J 2 = 0.000625) at T = 0.1 we are above the tran- 
sition and observe the Bow-tie structure, from there the 
pinch points develop when lowering this temperature to 
around T = 0,005. For larger \J 2 \ the temperature 



T = 0.1 is already within the AFM phase, therefore one 
sees the development of the pinch points already at this 
temperature. This means that aside of very low tem- 
peratures the additional J 2 term seem not to have much 
impact on the correlations. 



V. CONCLUSION 

To sum up, we have extended earlier work on the 
classical Heisenberg model on the Kagome lattice, by 
studying the inclusion of next-nearest neighbor interac- 
tions. This is relevant for real systems and a good way 
to gain further insight to the nearest neighbor model. 
We see that tiny values of J 2 have a big impact on this 
model by enabling finite-temperature phase transition, 
with larger transition temperature compared to the 
crossover point for J 2 = 0. This suggests that in real 
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materials having large spin on this lattice AFM order 
will be observed. The nature of the phase-transition 
has been explored, and it was found to fit the 2-D Ising 
universality-class for | J2I > 0.01, and shows consistency 
with a first order phase transition at the negative branch 
close to J 2 = 0. As the underlying induced symmetry 
is Z2 an Ising transistion appears reasonable. How 
this transistion changes towards first-order might be 
an inter sting point for future studies. For the tiny and 
positive J 2 case, two peaks in the specific heat develop, 
the second (in terms of higher temperature) most likely 
connecting further to the crossover point for the J2 = 
case. This is an interesting point for the nearest neighbor 
model and a possible line of attack for future studies. 
We notice that the lower transition line for the positive 



J2 branch between the AFM 2 state and the disordered 
phase tend to T — » for J 2 — » 0, while in the other 
branch the AFM 1 order seems to vanish before reaching 
the J 2 = point. This is consistent with earlier claims 
of a T = critical point for the AFM 2 transition. 
There seem to be no qualitative impact on the spin 
structure factor when tuning J ^. The region very close 
to J2 = remains still very interesting and should be 
further studied by analytic and numeric methods. 
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